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ABSTRACT 

Circumbinary disks are considered to exist in a wide variety of astrophysical objects, e.g., young 
binary stars, protoplanetary systems, and massive binary black hole systems in active galactic nuclei 
(AGNs). However, there is no definite evidence for the circumbinary disk except for some in a few 
young binary star systems. In this Letter, we study possible oscillation modes in circumbinary disks 
around eccentric and circular binaries. We find that progarde, nonaxisymmetric waves are induced in 
the inner part of the circumbinary disk by the tidal potential of the binary. Such waves would cause 
variabilities in emission line profiles from circumbinary disks. Because of prograde precession of the 
waves, the distance between each component of the binary and the inner edge of the circumbinary disk 
varies with the beat period between the precession period of the wave and the binary orbital period. 
As a result, light curves from the circumbinary disks are also expected to vary with the same period. 
The current study thus provides a new method to detect circumbinary disks in various astrophysical 
systems. 

Subject headings: black hole physics - accretion, accretion disks - binaries:general - stars: formation 
- planetary systems: protoplanetary disk - galaxies: nuclei 



1. INTRODUCTION 

Astrophysical disks are ubiquitous in the various sys- 
tems of the universe: star-compact object systems, star- 
planet systems, active galactic nuclei (AGNs), and so 
forth. These disks surround the individual objects as a 
circumobject disk. If the object surrounded by a rotat- 
ing disk is a binary, the disk is called a circumbinary disk 
(see Figure [T] for a schematic view of the circumbinary 
disk). 

About 60% of main sequence stars are consid- 
ered to be born as bin ary or multiple systems 
(|Duauennov fc Mavorl [T991). Numerical simulations 
have confirmed that young binary stars embedded 
in dense molecular g a s have a circumbinary disks 
(lArtvmowicz fc Lubowl Il996at iBate fe Bonnelil 119971 : 
iGimther fc Kiev! l2002t l2004f ) and that a circumsteller 
disk is also formed around each star (jBate fc Bonnelil 
[19971 IGimther fc Klevll200"l 1207)1 . Indeed, direct imag- 
ing of the circumbinary disk was successfully achieved by 
interferometric observat ions of a few youn g binary sys- 
tems including GG Tau (jDutrev et al.lll994[) and UY Aur 
(puvert et alj|1994f ). 

In the early stage of planet formation, a planet orbit- 
ing a star will be embedded in a rotatin g disk (hereafter 
circum binary disk) surrounding them. iKlev fc Dirksenl 
(2006) showed, by performing numerical simulations, 
that the cicrumbinary disk becomes eccentric due to the 
resonant interaction between the planet and the disk. 
Such a planet-disk interaction also causes the signif- 
icant evolution of the orbital elements of the planet, 
which gives a possible explanation about the observed 
high orbital eccentricitie s in extrasolar planetary systems 
(jGoldreich fc SarT1l2003( ). The direct probing of the cir- 
cumbinary disk is, therefore, a key to testing this sce- 
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nano. 

Massive black holes are c onsidered to co - evolve 
with their host g alaxie s (iMagorrian et all 119981 : 
iFerrarese fc Merrittl 120001: iGebhardt et all l2000ft . 
There is inevitably an evolutionary stage as a binary 
black hole in the course of gal axy merger until the 
coalescence of two black hol e s (iBegelman et all 119801 : 
iMaver et all 120071: lHavasakil 12008ft . lHavasaki et al l 



( 2007ft found that if a binary black hole is surrounded 
by a circumbinary disk, mass is transferred from the 
disk to each black hole. The system then has a triple 
disk composed of an accretion disk around each black 
hole and a c ircumbinary disk as a mass reservoir around 
the binary (|Havasaki et al.l [2008D . There is, however, 
still little evidence for the circumbinary disk as well as 
for the binary black hole itself. Therefore, the detection 
of the circumbinary disk is also one of the important 
scientific motivations in probing massive binary black 
holes with parsec/ s ubpar sec se parations. Q u ite rec ently, 
iBogdanovic et al.l (|2008ft and iDotti et all (|2008ft pro- 
posed the hypothesis that SDSSJO92712.65+294344.0 
is a massive binary black hole, by interpreting the 
observed emission line features as those arising from the 
mass-transfer stream from the circumbinary disk. 

There are many phenomena caused by oscillations in 
circumstellar/accretion disks. One of the most famous 
among them is superhumps. A superhump is a periodic 
luminosity hump on the light curves in an accreting bi- 
nary system, with a slightly longer period than the or- 
bital period of the binary. The superhump phenomenon 
was first discovered in the SU Ursae Majoris class of 
dwarf novae, which consists of a white dwarf and a late- 
type star with a low mass ratio, less than 0.2 (Patterson 
1979; Vogtl ll980f ). It is attributed to the precession of 
a deformed accretion disk i nduced by the tidal potential 
of the binary (|Osakilll985ft . iLubow I (|1991ft showed that 
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Primary Object Secondary Object 



Circumbinary Disk 

Fig. 1. — Schematic diagram of a circumbinary disk surrounding 
the primary object and the secondary object which are gravitation- 
ally bound as a binary. 

the deformation of the disk is due to the growth of an 
eccentric (i.e., m = 1) perturbation through nonlinear 
coupling with the tidal potential. Hitherto, no phenom- 
ena caused by oscillation modes have been detected in 
circumbinary disks. 

In this Letter, we investigate the tidally induced os- 
cillation modes in circumbinary disks. These modes can 
be driven by resonantly excited modes at particular reso- 
nance radii, which are similar to the ones responsible for 
superhumps in dwarf novae systems. The Letter is orga- 
nized as follows. In Section 2, we derive the azimuthally 
and temporally averaged tidal potential around an ec- 
centric binary and discuss the possible oscillation modes 
and their frequencies in circumbinary disks. Section 3 is 
devoted to a summary and discussion. 

2. TIDALLY INDUCED SPIRAL WAVES 

We consider two gravitationally bound point masses 
surrounded by a circumbinary disk. The point masses 
follow Kepler's third law, while the geometrically thin 
circumbinary disk has a rotation slightly deviated from 
Keplerian rotation because of the tidal potential of the 
binary. We assume that the circumbinary disk is trun- 
cated at an inner radi us, r- mi by the resonant intera ction 
with the binary (e.g., lArtvmowicz fc Lubowlll994l ). but 
extends outward to a large distance. 

2.1. Tidal potential 

First, we consider the tidal potential around an eccen- 
tric binary. In the inertial frame of reference centered 
on the common center of the binary, the gravitational 
potential can be written as 

H r M =- Gm GM 2 
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where M\ and M 2 are the masses of the primary and 
secondary components, respectively, ri and r 2 are their 
position vectors, and / is the true anomaly of the primary 
component. Here, {i — 1,2) periodically varies with 
time because of the orbital eccentricity and is given by 
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where a is the semimajor axis and 771 and n 2 are q/(l + q) 
and 1/(1 + q) with mass ratio q = M 2 /Mi, respectively. 



Azimuthally and temporally averaging equation {T]) 
with an assumption of r > ri and r 3> r 2 , we obtain 
the averaged binary potential over one orbital period, 
<t(r), in the form: 
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M 2 . Here, we have 
neglected terms of the order of e 3 or higher and made 
use of the following expansion formula: 



1 — e cos fl or ht H (1 — cos 2f2 or bt) + 0(e 3 ) 



(e.g.. iMurrav &: Dermottj 119991 ). Note that these series 
converge for e < 0.66. 

The second term on the right-hand side of equation §3§ 
gives the tidal perturbation potential around the eccen- 
tric binary: 
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From this equation, equation ([3]) is rewritten as $(r) — 
—GM/r + ^tidM- It i s noted from equation (0| that 
the tidal force acts more strongly on the circumbinary 
disk around an eccentric binary than around a circular 
binary. 

2.2. Trapped waves 

Let us now consider oscillations in the form of the nor- 
mal mode which varies as exp[i(ujt — k r r — m<f)\ with u> 
being the oscillation frequency and k T and m being the 
radial and azimuthal wave numbers, respectively. Then, 
the local dispersion relation for these oscillations is given 

by 
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where c s is the sound speed, O is the rotation frequency 
of the circumbinary disk, and k is the epicyclic frequency 
given by 



= 2ft(2f2 + rdn/dr). 
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From equation ([5]) , the region where a wave can prop- 
agate, i.e., k 2 > 0, is given by 
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In other words, waves can propagate inside the inner 
Lindblad resonance (ILR) radius, where ui/m = Q — K/m, 
or outside the outer Lindblad resonance (OLR) radius, 
where uj/m = fi + ft/m. However, we are interested only 
in the former type of waves with < u>/m < Q — n/m, 
because they are trapped in the region between the inner 
radius of the circumbinary disk and the ILR radius, and 
are thus possibly excited. All other waves are running 
w aves and wil l never be excited. 

iKatol (|1983f ) pointed out that the only possible global 
waves in nearly Keplerian disks such as circumbinary 
disks are low-frequency, one-armed (i.e., m = 1) waves. 
Waves with m ^ 1 have wavelengths much smaller than 
r, and are expected to damp rapidly as a result. In the 
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case of a circumbinary disk around a binary with nonex- 
treme mass ratio, however, mass transfer occurs from two 
points of the inner edge of the cir cumbinary disk to eac h 
component of every binary orbit (|Havasaki et aLlfeOOTl) . 
This might suggest that two-armed (i.e., to — 2) spi- 
ral waves are excited in the circumbinary disk through 
coupling with the tidal potential. 

2.3. Precession frequency and beat frequency 

From equation ([3]) , the orbital angular frequency of the 
circumbinary disk under the tidal perturbation potential 
is given by 
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where Ok = \/ GM/r 3 is the Keplerian angular fre- 



quency of the circumbinary disk, 
tion ((8]) in equation ((6]), we have 
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Substituting equa- 
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In general, the eigenfrequency and the extent of the prop- 
agation region can be estimated, using the WKBJ ap- 
proximation, by searching for oj that satisfies J k r dr ~ 
n 7r with an integer n, where the integration is performed 
over the propagation region. Given that the pressure 
gradient restoring force is much smaller than the grav- 
ity in geometrically-thin circumbinary disks and that the 
waves are trapped in the inner part of the disk, the pre- 
cession frequency of an m-armed spiral wave, ui Pim , is 
approximately given by 



(10) 



(11) 



at the inner disk radius, r = r; n . More rigorously, in the 
case of a one-armed (m = 1) wave, one can show that 
the precession frequency is significantly lower than, but 
is still of the order of, 0(ri n ) — K(r- m ). Thus we write it 

as 

< 3 ( a y /2 g d i y 

O orb ~ 4 \r iu ) (l + q) 2 \ 2 

Note that for m = 1 modes the precession frequency is 
positive and much lower than the orbital frequency. This 
feature means that the to — 1 waves with wavelengths 
comparable to r; n are trapped in the inner part of the 
circumbinary disk. 

On the other hand, the precession frequency of an m- 
armed (m ^ 1) wave normalized by the orbital frequency 
of the binary is obtained as 
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Thus, the precession frequency of m ^ 1 waves is approx- 
imately given by (m — l)QK.( r m)- Note that these waves 
have wavelengths much shorter than those of m = 1 
waves. 



It is instructive to evaluate the beat periods for some 
important modes. The beat frequency for an m-armed 
spiral wave is defined by 

Obcat.m = O orb — ^p,m- (13) 

From equation (|13p . we obtain the beat period as 
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According to the tid al truncation theory 
(|Artvmowicz fe Lubowl I1994D . the circumbinary disk 
around the binary with a low-moderate orbital eccen- 
tricity is truncated at r; n ~ 2.08a. In an equal-mass 
binary (q = 1) with e = 0.5, the beat period of the 
m = 1 mode is Pbeat.i ~ 1.03P or b, whereas that of the 
m = 2 mode is P bo at,2 ~ l-34P orb . 

3. SUMMARY AND DISCUSSION 

We have studied possible oscillation modes in a cir- 
cumbinary disk induced by the tidal potential of a bi- 
nary system such as young binary stars, protoplane- 
tary systems, and massive binary black holes in AGNs. 
We have pointed out that observationally interesting 
waves are those with precession frequency Lo p ^ m /m < 
O(nn) — K(r; n )/m, where r ln is the inner radius of the 
circumbinary disk. These waves are trapped between the 
inner disk radius and the ILR radius, where the pattern 
speed of the wave is equal to O(r) — n{r)/m. Among 
them, only m = 1 waves have wavelengths comparable 
to the inner disk radius. When the m = 1 mode is dom- 
inant, the inner region of the circumbin ary disk becomes 
eccentric and precesses very slowly (c.f., Kiev &: Dirksen I 



120061 iMacFadven fc Milosavlievicll2008h . On the other 
hand, waves with m ^ 1 have much shorter wavelengths 
and affect only the innermost narrow region of the disk. 
For example, if the m = 2 mode is dominant, the disk 
inner edge is deformed to an elliptical shape, which pre- 
cesses at the local Keplerian frequency. 

It is important to note that there are two types of 
excitation mechanism for these waves. In circular bina- 
ries, nonaxisymmetric perturbations in the disk can grow 
through the resonant interaction with the tidal potential 
at particular resonance radii, a mechanism similar to that 
for superhumps in dwarf novae systems. For example, 
the growth of an eccentric perturbation is driven by the 
excitation of an to = 2 wave at the 1:3 OLR radius. In 
addition to this mechanism, in eccentric binaries, a one- 
armed (m = 1) spiral wave can also be excited through 
dire ct driving as a result of a one-a rmed (m = 1) poten- 
tial (jArtvmowicz fc Lubow I[l996b| ). 

Whatever the excitation mechanism, the deformed in- 
ner part of the circumbinary disk precesses at the fre- 
quency given by equation (jTTJ) for the to = 1 mode and 
equation (fT2|) for the m = 2 mode. Since the veloc- 
ity field in the disk is also perturbed by the waves, the 
emission line profiles from the inner part of the circumbi- 
nary disk will vary with the precession period. Such a 
variability has a distinct feature and will easily be ob- 
served. Since the m = 1 mode is a low-frequency, ec- 
centric mode, the relative intensity of the violet (V) and 
red (R) peaks of double-peaked line profiles varies with 
a long period, e.g., ~ 40P or b for an equal-mass binary 
with e = 0.5. Such a line-profile variability caused by 
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m = 1 waves has long been known as the V/R varia- 
tion in Be stars , B-type stars with circ umstellar decre- 
tion disks (e.g.. iPorter &: Riviniusl [20031 ). In contrast, if 
the m = 2 mode is dominant, the double-peaked pro- 
files stay symmetric, but their peak separations (and 
FWHMs) will vary with a short period of ~ 2 3 / 2 P orb 
irrespective of binary parameters. 

In addition to these line-profile variabilities, another 
type of variability from circumbinary disks is expected. 
There is the beat between the orbital frequency of the 
binary and the precession frequency of the circumbinary 
disk. The beat period is slightly longer than the orbital 
period for m = 1 and m = 2 modes, as seen from equa- 
tions (jXSJ) , and (jT4j) . The radiation emitted from 
the circumbinary disk through the tidal dissipation is 
expected to vary periodically with the beat period, be- 
cause the distance between each component of the bi- 
nary and the inner edge of the circumbinary disk period- 
ically changes. Moreover, the mass transfer rate from the 
circumbinary disk to each binary component is also ex- 
pected to vary with the same period, because of angular 
momentum removal by the tidal torques. Therefore, the 
circumbinary disk perturbed by a trapped density wave 
will show the light-curve modulation at the beat period. 

As for observability, the light-curve modulation caused 
by an m = 1 deformation is expected to be detected more 
easily than those caused by the m ^ 1 one. This is bc- 
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cause the m = 1 wave can exist globally in the circumbi- 
nary disk, while the m ^ 1 waves exist only in a narrow 
region. 

In this Letter, we have shown that variabilities with the 
beat period and/or the precession period are expected as 
a natural consequence of the tidal interaction between a 
binary and a circumbinary disk and that one can identify 
circumbinary disks in terms of these periodic variabili- 
ties. This provides a new method to probe circumbinary 
disks in various astrophysical systems. In a forthcoming 
paper, we will perform a more detailed analysis of the 
mode characteristics, including numerical simulations. 
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